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We study the electron-energy loss spectra of strongly correlated electronic systems doped away 
from half-filling using dynamical mean-field theory (d = oo). The formalism can be used to study 
the loss spectra in the optical (q = 0) limit, where it is simply related to the optical response, and 
hence can be computed in an approximation- free way in d = oo. We apply the general formalism 
to the one-band Hubbard model off n = 1, with inclusion of site-diagonal randomness to simulate 
effects of doping. The interplay between the coherence induced plasmon feature and the incoherence- 
induced high energy continuum is explained in terms of the evolution in the local spectral density 
upon hole doping. Inclusion of static disorder is shown to result in qualitative changes in the low- 
energy features, in particular, to the overdamping of the plasmon feature, resulting in a completely 
incoherent response. The calculated EELS lineshapes are compared to experimentally observed 
EELS spectra for the normal state of the high- Tc materials near optimal doping and good qualitative 
agreement is found. 



PACS numbers: 75.30.Mb, 74.80.-g, 71.55.Jv 

The ac conductivity and dielectric tensor provides 
valuable information concerning the finite frequency, fi- 
nite temperature charge dynamics of an electronic fluid 
in a metal. The dramatic changes in the electronic state 
from localized to itinerant across an insulator-metal (I- 
M) transition, be it driven by pressure, or doping the in- 
sulator, is reflected in concomitant changes in the optical 
responses. These studies, therefore, give one a systematic 
picture describing the nature of the change in electronic 
state, appearance of new low-energy excitations, their 
dispersion and stability. In cuprate superconductors, for 
e.g, these studies have convincingly demonstrated the 
non-Fermi liquid nature of the charge [Q and spin 
dynamics in the normal state. 

Optical response in a solid provides an estimate of 
the carrier scattering rate at finite frequency, and thus 
give detailed information about the finite frequency col- 
lective excitations responsible for scattering the carriers. 
In a quantum paramagnetic metal, the dominant carrier- 
scattering mechanism are electron-hole pairs, which are 
the low-energy collective excitations- it is precisely the 
e-h spectrum which is measured by electron-energy-loss 
spectroscopy (EELS), which provides one with a direct 
experimental picture for the spectral density of particle- 
hole excitations in a solid 

In a weakly interacting system, one expects the details 
of the particle-hole spectrum to be sensitive to details of 
the shape and size of the Fermi surface (e.g, its curva- 
ture). That such a picture is untenable for strongly cor- 
related systems was pointed out by Shastry ei a/ Q , who 
argued that in this case, the whole Brillouin zone tends 
to get populated. Another important point of consider- 
able relevance is the transfer of spectral weight over large 
energy scales that is characteristic of strongly correlated 
systems-it has no analog in weakly interacting fermi sys- 



tems. Consequently, one expects qualitatively different 
responses in this case, compared to those characteristic 
of weakly interacting systems. 

Quite generally, the transmition EELS lineshape is re- 
lated to the wave- vector and frequency dependent dielec- 
tric function via the equation, 

/(q,.)^-/..-^. (1) 

Following Rcf. we employ the random-phase- 

approximation (RPA) to treat the small-q part. In 
the optical limit that we are interested in {q 0), 



e(q ^ 0,t^) = l-fdTT 



where (t(u;) is the longitudi- 



nal optical conductivity. The analysis carried out along 
these lines is invalid when ga ~ 1, where the short-ranged 
part of the potential has to be considered and the RPA 
is inadequate. However, local correlation effects are cor- 
rectly incorporated into the above eqn via cr(u;), which 
can be computed reliably within d = oo, for example. 
In this limit, q — 0, the particle-hole spectral density is 
obtainable from the dissipative part of H , 



1(0, Lu) = -Im 



1 



(2) 



with e'(w) = 1 - ^^^^ and e" (w) = go that one 

needs to have an estimate of the longitudinal dielectric 
constant to calculate the EELS spectrum. It is to be 
noted that this yields the EELS lineshape corresponding 
to the experiment performed in the so-called "transmis- 
sion mode" . Thus, the problem has now been reduced 
to that of computing e(a;), which, being a local quan- 
tity, can be computed in an approximation-free way in 
d = oo ||. 

In this communication, we want to develop a theory 
of the trasmition EELS (eq. (|l|)) for electronic systems 
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which undergo Mott transitions as a function of band- 
filling 1^. To be specific, we consider the one-band Hub- 
bard model, where the Mott transition can be driven by 
hole-doping the Mott insulator, or by pressure (which 
changes the ratio U/t). To capture more fully the effects 
of chemical substitution- (which is how doping is carried 
out in practice), we also introduce static, site-diagonal 
disorder in the model. The resulting hamiltonian is, 

iJ = - ^ tijc\^c.ja + ni^nii - ^(wi - p)nia (3) 



as a prototype model describing the electronic degrees of 
freedom in TM oxides. To study the 3D case, we em- 
ploy the d — oo approximation, which is the best among 
those currently available to study the M-I transition 
Since this method has been extensively reviewed, we only 
summarize the relevant aspects. All transport properties, 
which follow from the conductivity tensor, are obtained 
from a k-independent self-energy in d = c>o; the only in- 
formation about the lattice structure comes from the free 
band dispersion in the full Green fn: 



UJ + l^i - Ek - S(w) 



(4) 



To solve the model in c? = oo requires a reliable way to 
solve the single impurity Anderson model(SIAM) embed- 
ded in a dynamical bath described by the hybridization 
fn. A(ti;). There is an additional condition that com- 
pletes the self consistency: 



deG{e,iuj)pa{e) = 



1 



(5) 



iLu + ji — A{iuj) — S(iw) 

where po(e) is the free DOS {U — 0). The above eqns. (^ 
^) refer to the disorder-free Hubbard model. With micro- 
scopic, binary disorder (which is the only type we con- 
sider here). 



P{v,) = (1 - x)S{v,) + x6ivi - v) 



(6) 



only the disorder-averaged quantities are physically ob- 
servable. In this case, we use an extended version of 
the iterated perturbation theory (IPT), which combines 
the effects of dynamical correlations (via usual IPT) with 
those of disorder (via CPA) in a selfconsistent way. For 
the sake of completeness, we outline the calculational 
details briefly. As a first step, we compute the full lo- 
cal Green function for the pure model using IPT. This 
has been shown to yield results in excellent agree- 
ment with those obtained from exact diagonalization 
techniques. This IPT propagator is then corrected for 
repeated scattering from local disorder by computing the 
new self energy from the usual CPA §, 



-(1 - x)Y.{uj) 
1 + Y.{lu)G{lo) 



x{v — S(ij-')) 
1 - (i; - S(w))G(cj) " 

(7) 



which results in an implicit equation for the interaction 
and disorder corrected self energy (and Green function). 
To make the treatment selfconsistent, we compute the 
new IPT selfenergy using the new GF computed with 
(IPT+CPA) in the first step. This procedure is repeated 
until convergence is achieved. At each step of the iter- 
ation, we fix the chemical potential from the Luttinger 
sum rule, n = J!^^ p{uj)duj = x. This extended IPT 
yields the local propagator, < Gii{uj) > corrected both 
for dynamical correlations and disorder induced repeated 
scattering, both treated on the same footing In 
d — oo, this is sufficient to compute the transport, be- 
cause the vertex corrections in the Bethe Salpeter eqn. 
for the conductivity vanish identically in this limit [Q. 
Thus, the conductivity is fully determined by the basic 
bubble diagram made up of fully interating local GFs of 
the lattice model. 

The optical conductivity and the Hall conductivity are 
computable in terms of the full d = oo GFs as follows : 

o-xx{iu!) =^ — depo{e)y^G{e,iiy)G{e,ii^ + iuj) (8) 

lUJ J ^ 

The dynamical dielectric constant is directly calculated 
from the optical conductivity via. 



47r 

^xx{^) = H i'7^2;(w) 



(9) 



yielding its real and imaginary parts. Use of eqn. (3) then 
permits us to study Ieels{^) ^^^7^^ 

Before embarking on our results and their analysis, it 
is instructive to summarize what is known about the 
d — oo Hubbard model. At large U/t, and away from 
half-filling (n = 1), the ground state is a paramagnetic 
FL if one ignores the possibility of symmetry breaking 
towards antiferromagnetism, as well as disorder effects, 
which are especially important near n — I. This can be 
achieved formally by introducing a nnn hopping, which 
in d = oo leaves the free DOS essentially unchanged 
This metallic state is characterized by two energy scales: 
a low energy coherence scale Tcoh, below which local FL 
behavior sets in and a scale of 0{D), {D is the free 
bandwidth) characterizing high energy, incoherent pro- 
cesses across the remnant of the Mott-Hubbard insulator 
at n = 1. At r < Tcoh, the quenching of the local mo- 
ments leads to a response characteristic of a FL at small 
Lo « 2D (but with the dynamical spectral weight trans- 
fer with doping, a feature of correlations), but at higher 
T > Tcoh, the picture is that of carriers scattered off by 
effectively local moments, which makes the system essen- 
tially like a non-FL (note that the metal with disordered 
local moments is not a FL). 

The picture is modified drastically in presence of static, 
diagonal disorder, in which case, the extended IPT has 
to be used. For weak disorder, the above picture is qual- 
itatively unaffected, and the only interesting features are 
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that Luttinger's theorem is violated, and the FL quasi- 
particles acquire a finite Hfetime at [^. However, with 
increasing disorder, repeated scattering effects destroy 
the quasiparticle picture, and the metalhc state is charac- 
terized by an incoherent, non-FL response with a pseudo- 
gap structure in the DOS. Similar results have been ob- 
tained earlier in Ref . Q , where a T- matrix approximation 
was employed to treat disorder-induced resonant scat- 
tering. The above describes the competing tendency of 
the correlation-induced low temperature coherence scale, 
Tcoh in the d = oo Hubbard model, with the disorder- 
induced incoherence, which tends to suppress it, driving 
it eventually to zero, leading to a non-FL metal state. 
A large enough disorder strength leads to band-splitting, 
as in usual CPA, leading to a continuous transition to a 
disorder driven insulating state. 
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FIG. 1. Real part of the optical conductivity (t^^ for 
U/D = 3.0, (5 = 0.1 (continuous) and S = 0.2 (dot-dashed). 
The lower panel shows the imaginary dielectric constant for 
the same parameter values. Note the isosbectic point in both 
quantities. 

Armed with this information, we are ready to discuss 
our results. We choose a gaussian unperturbed DOS, and 
U/D = 3.0 to access the strongly correlated FL metal- 
lic state off n = 1, ignoring the possible instability to 
an AF-ordered phase. All calculations are performed at 
a low temperature, T = O.OII?. We are mainly inter- 
ested in the variation of the EELS lineshape with hole 
doping, given here hy S ~ {1 — n). This fixes the chem- 
ical potential, and the FL resonance position, and the 
IPT describes the evolution of spectral features in good 
agreement with exact diagonalization studies Q . In view 
of the ability of the IPT to reproduce all the qualitative 
aspects observed in (Txx{i^), we believe that is a good 
tool in the present case. Fig.(l) shows the optical con- 
ductivity and the longitudinal dielectric constant, a^xi^) 
agrees with calculations performed earlier in all the main 



respects; in particular, it clearly exhibits the low-energy 
quasicoherent Drude form, the transfer of optical spectral 
weight from the high-energy, upper-Hubbard band states 
to the low energy, band-like states with increasing hole 
doping, and the isosbectic point at which all the a^xii-^) 
curves as a fn. of filling cross at one point, to within 
numerical accuracy. It is interesting to point out that 
such features have also been observed in experimental 
studies Correspondingly, Imexxi^^) ^-Iso shows the 

isosbectic point, the explanation for which is identical to 
the heuristic one proposed recently by us for the case of 
axx{uj) 0. 




FIG. 2. The electron-energy loss spectra for the doped 
Hubbard model for U / D = 3.0, and 5 = 0.1, 0.2, 0.3 as shown. 
The curve with circles denotes the Drude contribution from 
renormalized single-hole excitations in the metallic phase. 

In Figs. 1^ and |^, we show the theoretically calculated 
EELS spectrum for U/D = 3.0 and for different band- 
fillings as a function of disorder strength, v. This al- 
lows us to access the intermediate correlation regime, 
where the absence of any small parameter in the prob- 
lem requires use of controlled methods that interpolate 
correctly between the weak- and strong coupling limits 
(such as the d = oo method used here). From the cal- 
culation of the optical response, we have checked that 
luYUxxi^^ = 0) = 0, implying that — i?el/e(0) = 0, as 
required ||l2|. With U/D — 3.0, we observe a peak at 
Lo ~ U /2, associated with the transitions from the lower- 
Hubbard band to the quasiparticle resonance. We as- 
sociate this peak to the strongly renormalized particle- 
hole excitations in the strongly correlated metal. In this 
context, two features are worthy of mention: at low en- 
ergy, the loss intensity goes quadratically with uj, and the 
particle-hole peak position shifts with increasing doping. 
The shift of this feature correlates with that of the mid- 
IR peak in the optical conductivity (Fig.(l)), and is thus 
a clear manifestation of the transfer of optical weight 
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to lower energy with hole-doping. Moreover, appreciable 
loss intensity is observed at higher energies (w > 2.0) as a 
smooth continuum, in contrast to the expectations from 
a weakly correlated (very small U/D) FL. This clearly 
shows how the delicate balance between coherence and 
incoherence is controlled by the transfer of optical spec- 
tral weight in the full EELS spectrum. In fact, we ex- 
pect the coherent feature in the full spectrum to become 
sharper as T (O.OID in our analysis) is lowered further. 
We see from the above that coupling to high-energy, in- 
coherent, multiparticle excitations strongly renormalizes, 
but does not destroy, the coherent response in a strongly 
correlated Fermi liquid. However, the analysis presented 
above shows that a proper treatment of coherence (re- 
lated to itinerance) and incoherence (coming from the 
local, atomic-like features) on an equal footing is neces- 
sary to obtain consistent results. The EELS at very low 
energy are well described by a Drude fit, but features like 
the damped plasmon peak at uj/D ~ 1.0 are distinctly 
non-Drude-like, and can only be reliably accessed by the 
full analysis, as described above. 
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FIG. 3. The EELS spectra for the d = cx) Hubbard model 
for U/D = 3.0 and v/U = 1/4, 1/3, 1/2. See the text for the 
explanation of various structures in I eels (i^) in the light of 
the underlying structure of the local spectral density of the 
HM in d = OQ. 

Given the above, we expect small coherence destroying 
perturbations to tilt the balance in favor of low-energy 
incoherence. Consideration of doping-induced static dis- 
order has precisely such an effect on the low energy EELS 
spectrum. Consistent with the disorder-induced smear- 
ing of the Drude part in axxi^o), the coherent part in the 
EELS spectrum is replaced by a broad peak, which still 
shows remnants of the coherent peak at f = (u = J7/4) . 
However, the low-energy dependence in the pure case 
is replaced by a linear lu dependence in this case. The 
complete destruction of the quasiparticle behavior in the 



DOS and the Drude-like response at t; = J7/3 is con- 
sistent with the complete destruction of the "coherent" 
peak in the EELS spectrum; we understand this as over- 
damping of the collective p-h peak by disorder-induced 
strong scattering as v is increased. In the strong disorder 
regime, the metallic state is incoherent with a pseudogap 
feature in the DOS and the action of the potential 
Vkk' between the Hubbard band states does not create 
well-defined elementary excitations, leading to an inco- 
herent response at low energy. 

In experiments carried out on the cuprates in their 
normal state, a very broad plasmon peak is revealed [ p^ , 
implying poorly defined plasmon excitations in the cor- 
related metallic state. These studies also reveal that the 
optic "plasmon" disperses quadratically (in q). Recent 
studies indicate an acoustic-like heavily damped plas- 
mon mode at small q ||l3| with a linear oj dependence 
in Bi2Sr2CaCu20s, alongwith a single broad peak. We 
cannot directly compare our results to those of , the 
EELS experiment was carried out in the reflection mode, 
and this measures the quantity Img{Q,uj), which is not 
directly related to what we have calculated here. How- 
ever, looking at the results of we see that the es- 
sential features seen there are indeed reproduced by our 
calculation. To make contact between our calculations 
and experiments, we shall suppose that the one-band 
Hubbard model is understood to be the effective model 
describing the coupled spin-carrier dynamics in the Cu- 
O planes. We will also assume that bandstructure ef- 
fects are not crucial to the case under study, and that 
the source of the anomalous behaviors is rooted in the 
dynamical effects of the strong, local Hubbard interac- 
tion, an assumption well-supported at doping levels near 
the optimal We stress that the above conditions 
will be modified in the pseudogap phase, where the pre- 
cursor effects of d-wave pairing fluctuations-intrinsically 
non-local, are appreciable |15| ]. In this doping range, 
the dynamical effects of non-local correlations will prob- 
ably lead to even stronger deviations from a FL picture, 
though a concrete calculation including such effects re- 
mains to be developed. With these caveats in mind, 
we find good qualitative agreement between our results 
for V = t/73 of Fig. ^ and the experimental EELS line- 
shapes [O]; the quadratic (in lu) behavior at low w, the 
broad "plasmon" peak (strongly damped) , the asymmet- 
ric lineshape and the continuum response at higher en- 
ergy are all reproduced in qualitative agreement with ex- 
periment. It would be interesting to check whether the 
plasmon peak shifts to lower energies with hole doping; 
this would be an interesting check on the validity of the 
approach presented here. The comparison is quite good 
upto OJ / D ~ 3.0, beyond which the simple single-band 
modelling does not apply anyway. 

There is extensive experimental work [ p^ indicating 
that the metallic phase above near optimal doping in 
cuprates is not describable in terms of Landau Fermi liq- 
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uid ideas. Prelovsek et al. have recently considered 
the question of the EELS hneshape in the t — J model, 
making use of results obtained from finite-temperature 
Lanczos techniques. Within our approach, it is known 
that the metallic state off n = 1 above the lattice co- 
herence scale, Tcoh, is not a Fermi liquid the cal- 
culation carried out above is then consistent with the 
non-FL charge dynamics if the Tcoh can be driven suffi- 
ciently low (below 0.0 ID used here). We expect precisely 
a very low Tcoh near n = I for U/D near the Mott transi- 
tion ^ . Comparing our results with those found in [ pT| , 
we observe that our results are in nice agreement with 
theirs, but the results presented here are more transpar- 
ent, being based on an analytical scheme. In addition, 
the results show that the d = oo approach is capable of 
producing the low- and high-energy features on the same 
footing, and moreover, permits us to include effects of 
doping-induced static disorder in a consistent way. In 
view of this, we believe that the approach described here 
should also be applicable to a wide variety of doped tran- 
sition metal oxides, where the combined effects of corre- 
lations and disorder are especially pronounced Ad- 
ditionally, effects of orbital degeneracy, etc, in real ma- 
terials can be treated by a suitable generalization of the 
above procedure This is a non-trivial problem, be- 
cause the solution of the impurity model is harder. We 
hope to address this issue in future work. 

In conclusion, we have considered the energy loss func- 
tion of a model for transition metal oxides, and captured 
the effects of strong, local correlations and static disor- 
der in a consistent way. Comparison of our results with 
the experimental EELS spectra for cuprates in the "nor- 
mal" state above Tc near optimal doping shows that all 
the main observed features are consistently reproduced 
by our calculation. 
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